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We investigate the particle production of a scalar field \ coupled to an inflaton field (f> (g 2 (f> 2 X 2 /2) 
in the oscillating inflation model, which was recently proposed by Damour and Mukhanov. Although 
the fluctuation of the <j> field can be effectively enhanced during a stage of the oscillating inflation, 
the maximum fluctuation is suppressed as the critical value <f> c which indicates the scale of the core 
part of the inflaton potential decreases, in taking into account the back reaction effect of created 
particles. As for the x particle production, we find that larger values of the coupling constant g are 
required to lead to an efficient parametric resonance with the decrease of 4> c , because an effective 
mass of inflaton around the minimum of its potential becomes larger. However, it is possible to 

■ generate the superheavy x particle whose mass is greater than f0 14 GeV, which would result in an 

' important consequence for the GUT baryogenesis. 
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I. INTRODUCTION 



Inflationary cosmology is one of the most reliable concepts to describe the early stage of the universe This 
i paradigm not only gives the solution to a number of shortcomings of the standard big bang cosmology, but also 
provides the density perturbation that may be responsible for the structure formation During the inflationary 
stage, a scalar field <p known as inflaton is slowly rolling down toward a minimum of its potential. Inflation ends when 
O ; the kinetic energy of the inflaton becomes comparable ,„ the potential energy. After that, the inflaton fleld begins to 
oscillate around the minimum of its potential and produces elementary particles. This process by which the energy of 
Oh, the inflaton field is transferred to other particles is called reheating. The original scenario of reheating was considered 

D 

1 field. However, since the production of a grand unified theory (GUT) scale boson is kinematically forbidden in this 



in Ref. H based on the perturbation theory, adding the phenomenological decay term to the equation of the inflaton 



model, the GUT scale baryogenesis does not work well in this scenario. 

It was recently recognized that the reheating process begins by an explosive particle production called preheating 
[Q,^). During this stage, the fluctuation of produced particles grows quasi-exponentially by parametric resonance. The 
efficiency of resonance depends on the model of the inflation. The chaotic inflation is one of the most efficient model 
for the development of the fluctuation. In the case of the massive inflaton potential V(4>) = m 2 <fr 2 /2, another scalar 
field x coupled to inflaton with a coupling g 2 (f) 2 x 2 /2 can be enhanced in a certain range of the coupling constant 
||-^]]. With the coupling of g > 10~ 4 , resonance turns on from a broad resonance regime and the fluctuation of the 
X particle increases overcoming the diluting effect by the cosmic expansion. In the massless inflaton potential X4> /4, 
the growth of the inflaton fluctuation occurs even if we do not introduce another field x coupled to inflaton J§-[t| . In 
this case, however, since the resonance band is restricted to be narrow, the transfer of energy from the homogeneous 
inflaton to the fluctuation does not occur sufficiently. In the two-field model of V{(f>, x) — A0 4 /4 + g 2 2 x 2 /2, the 
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X particle production is typically more efficient than the production of the (j> particle fWj. As for other models of 
inflation, several authors considered preheating such as the hybrid inflation jl^ , and the higher curvature inflation 



Recently it has been recognized by Damour and Mukhanov |17| that inflation occurs even in the oscillating stage 
for the non-convex type potential U(</>) where cPV/dcj) 2 is negative in the regime not too far from the core of the 
potential. In spite of the rapid oscillation of the inflaton field, inflation takes place during which the field moves in flat 
regions outside the core part. Liddle and Mazumdar [ fL8| called it oscillating inflation and numerically calculated the 
revised number of e- foldings. The amount of inflation becomes larger with the decrease of the critical value <j) c which 
determines the scale of the convex core of the potential. However, it is difficult to induce a sufficient inflation only 
by the oscillating inflation even when <j) c is lowered to the electro- weak scale (f> c = 10~ 17 Af p i p8|| . Since the ordinary 
inflation takes place in the slow-roll regime before the oscillating inflation, the total amount of e-folding is mostly due 
to this preceding inflation and the contribution of the oscillating inflation is added to some extent. The energy scale 
of the potential can be determined in the usual manner by fitting the density perturbation observed by the Cosmic 
Background Explorer (COBE) satellite. 

As was pointed out by Damour and Mukhanov, since the square of the effective mass of the inflaton field m 2 ^ = 
d 2 V/d<p 2 is mostly negative during the oscillating inflation, the fluctuation of inflaton can be strongly amplified. Taruya 
p9| considered the evolution of the fluctuation including the metric perturbation in the single field model during the 
oscillating inflation. Although the super-horizon modes (fc — > 0) are not relevantly amplified, the growth of the modes 
inside the Hubble horizon occurs significantly. In the case of the single field, since the curvature perturbation on 
the comoving slice remains constant in the long wavelength limit, we can not expect the significant growth of the 
long-wave perturbation. It was also suggested by Damour and Mukhanov that another field \ coupled to inflaton 
would be resonantly amplified during the oscillating inflation and the superheavy GUT scale bosons are expected to 
be generated. The structure of resonance for the x particle is different from that of the <j) particle, and we should 
make clear the parameter range of the coupling constant g where the x particle production occurs sufficiently. Since 
the effective mass of the inflaton around the minimum of the potential depends on the critical value 4> c , it is also 
important to investigate how the difference of the shape of the potential would affect the x particle production. If 
the x particle is efficiently created in the oscillating inflation model, this would affect a nonthermal phase transition 
ppf , and topological defect formation pi)) . Also, the possibility of the production of the superheavy particle would 
make the baryogenesis at the GUT scale possible p2| . In this paper, we consider the production of the x particle as 
well as the <f> particle during the oscillating inflation and the subsequent oscillating stage. We make use of the Hartree 
approximation in order to include the back reaction effect of created particles which is absent in Ref . |l9f| . Even in 
the case where the resonance band is very broad, the back reaction effect plays a crucial role to terminate the growth 
of the fluctuation. We will discuss how the maximum fluctuations of <fi and x particles depend on the critical value of 
4> c and the coupling constant g. 

This paper is organized as follows. In the next section, we introduce basic equations based on the Hartree approx- 
imation in the oscillating inflation model. In Sec. Ill, the particle production in this model is investigated by making 
use of the numerical results. We study how the fluctuations of the <f> and x fields grow during the oscillating inflation 
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and the subsequent oscillating stage. We give our discussions and conclusions in the final section. 



II. THE BASIC EQUATIONS 



We investigate a model where an inflaton field <f> is coupled to a scalar field \ p3fl , 
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(2.1) 



where k 2 /8n = G = M~ 2 is Newton's gravitational constant, R is a scalar curvature, m x is a mass of the x field, and 
<7 is a coupling constant. We consider an inflaton potential V(<p) proposed by Damour and Mukhanov ||I7|, which is 
described by 
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where M is a mass which is constrained by the primordial density perturbation observed by the COBE satellite, q is 
a dimensionless parameter greater than zero, and <p c is a critical value of the inflaton field which determines the scale 



of the core part of the potential (See Fig. 1). Note that the potential (2.2) is a toy model which is not directly related 
with a real theory of physics, although supergravity and superstring models may give rise to some non-convex type 



potentials jL7|. In the case of <fi 2 ^S> <fi 2 , the potential is approximately written by V(<p) ~ Aq 



) q , and d 2 V/d4> 2 



becomes negative for q < 1. Inflation takes place in the usual manner while inflaton slowly moves in this flat region. 
This conventional slow-roll inflation is followed by the oscillating inflation, which means that inflation continues to 
occur during the oscillating stage of inflaton while it evolves in the region of <fi 2 > <f> 2 . Since the amplitude of the 
<t> field gradually decreases due to the adiabatic expansion of the universe, the <f> field is finally trapped in the core 
region of the potential (i.e. cf> 2 < <fi 2 ) and the oscillating inflation ceases. After that, the universe enters the ordinary 



reheating stage, in which the potential (2.2) is described by the massive inflaton potential. 

Hereafter, we study the dynamics of the system in the flat Friedmann-Robertson- Walker metric 



ds 2 



-dt 2 +a 2 (t)dx. 2 , 



(2.3) 



where a(t) is the scale factor, and t is the cosmic time coordinate. 

Let us consider the equations of motion based on the Hartree factorization in the oscillating inflation model. We 
decompose the inflaton field into the homogeneous and fluctuational parts as 



0(t,x) = 0o(i) + $0(t,x), 
where the fluctuational part satisfies the tadpole condition 

(<ty(f,x)) = 0. 

In order to study the quantum particle creation, we expand 5<f) and \ fields as 
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X = 
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a>kXk{t)t 



where a k and a\ are the annihilation and creation operators respectively. 

Then, the equations of motion for inflaton are expressed by imposing the Hartree factorization 



as 
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where H = a/a, V n (<p ) = 5 n V((f>o)/6<f>%, and (<5</> 2 ), ( X 2 ) are defined by 

W> 2 ) = ^2 / k 2 \5fa\ 2 dk, 

(X 2 ) = ^J k 2 \ Xk \ 2 dk. 



(2.9) 
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Note that the V^ 2 \<Pq) term in Eq. (2.9) is the leading term for the development of the fluctuation (<5</> 2 ). In the 
present model, this term is negative in most stages of the oscillating inflation, and the growth of the fluctuation can 
be expected during this stage. As <fi particles are produced, however, the back reaction effect by the growth of (5<fr 2 ) 
plays an important role. The development of the fluctuation is suppressed by the (f> particle production itself. 
With regard to the x held, this satisfies the following equation 



Xk + "iHxk 
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Xk =0. 



(2.12) 



The g 2 (j)Q term leads to the parametric amplification of x particles when the fa field oscillates around the minimum 
of its potential. The strength of resonance depends on the coupling constant g and the amplitude of the fa field. 
Moreover, the effective mass m<f, of the fa field is also important as we will show later. Since is closely related 
to the shape around the core region, the development of the fluctuation (x 2 ) depends on the critical value (j> c . In 
the case where the growth rate of the x fluctuation is large, back reaction effects of x particles as well as <fi particles 
finally shut off the parametric resonance. 

Next, the evolution of the scale factor is described by 
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((Vx) 2 ) = / Axk?dk. (2.16) 
Before analyzing the above equation of motion, we first discuss the evolution of the scale factor and the homogeneous 



inflaton field neglecting the fluctuational terms. Then, Eqs. ( |2.q ) and (2.13) are approximately written as 

4>o + 3H4> + V,<f, (0o ) » 0, 

~0 O 2 + V{(f> ) 
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(2.17) 
(2.18) 



We can understand the mean behavior of the scale factor and the inflaton field by ignoring the core region of the 
potential. For the case of 0q 3> ^, making use of the time averaged relation (4> 2 )t = q(V(4>))T in Eqs. ( 2.17} ) and 



( 2.18 ), we find the following approximate relation |17J 
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where cf>o is the amplitude of the <fio field. Note that the universe expands acceleratedly for < q < 1. Hereafter, we 
consider the case of < q < 1. Since inflation takes place for 0q > (j) 2 , the amount of inflation becomes larger with 



the decrease of • 



Damour and Mukhanov estimated the number of e-foldings as 
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where the subscripts s and / denote the end of the slow-roll and the end of the oscillating inflation respectively. The 
slow-roll inflation ends when the slow-roll parameter e = (V,0 O /V) 2 /2k 2 becomes of order unity. For the case of 
0o ^ 0c) since e can be written as 
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4>o(t s ) is estimated by setting e = 1 as 
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We find from Eqs. (2.21) and (2.23) that the number of e-foldings during the oscillating inflation is small compared 
with the total needed amount of inflation N t > 60. For example, when cj) c = 10 _6 M p i and q = 0.1, N — 3.3. Even 
if <fi c is the electro- weak scale cj> c — 10~ 17 -M p i, N = 12.2 for q = 0.1. As was confirmed by numerical calculations in 
Ref. fl8[| , the number of e-foldings becomes smaller as q approaches zero or unity for the fixed value of <\> c . Hence 
even if we change the values of q in the range of < q < 1 , the contribution to the amount of inflation due to the 
oscillating inflation is still small. This means that the preceding inflation by the ordinary slow-roll is expected to 
contribute to most of the number of e-foldings needed to solve cosmological puzzles. 

Let us estimate the energy scale of the potential ( |2.2| ) which is constrained by the density perturbation observed 
by the COBE satellite. First, we consider the value of <j>Q (— <fio{ti)) at the epoch of horizon exit when physical scales 
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crossed outside the Hubble radius 50 e-foldings before the start of the oscillating inflation. Calculating the number 
of e-foldings 
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in the present model with the condition of ^ > ^, we obtain 
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Then the value of (f>o{ti) is approximately estimated as 
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where we neglected the contribution of the <f)o(t s ) term. The square of the amplitude of the density perturbation can 
be calculated for 3> <j> 2 as p5) 
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Then the mass M is constrained to be 
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The COBE data requires k 2 x 10 5 . Setting N = 50, we can estimate the mass M as two functions of q and 
(j) c . The mass M is weakly dependent on the value of <j> c for the fixed value of q. For example, for <p c — 10~ 4 M p i and 
q = 0.1, M = 3.74 x 10~ 4 M pl ; and for <\> c = 10" 6 Af pl and q = 0.1, M = 3.33 x 10~ 4 M pl . In the next section, we 



use the value of M which is obtained by Eq. ( 2.28 ), and consider the growth of the fluctuation during and after the 
oscillating inflation. 



III. PARTICLE PRODUCTION IN THE OSCILLATING INFLATION MODEL 



In this section, we investigate the particle production of 4> and \ fields with the potential (2.2). Defining new scalar 
fields (f = a 3 / 2 o , S(p k = a 3 ^ 2 S<pk, and X k = a 3 ^ 2 Xk, Eqs. (2.8), (2.9), and ( 2.12| ) can be rewritten as 
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where we have considered the 



terms up to n = 1 . The higher order terms do not significantly contribute to the 



evolution of the system. We find from Eq. (3.4) that the frequency of the tpo field changes even during one oscillation 
of the ipo field. At the stage when the particles are not sufficiently produced ((<5</> 2 ), (x 2 ) *C <Pq), the second and third 
terms in Eq. (3.4) are negligible. Moreover, since the last term in Eq. ( |3.4| ) is estimated by using Eqs. (2.17) and 
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this term is also negligible compared with the first term in Eq. (3.4) for the typical values of <j) c ( < 10 M p \). Then 



the effective mass of the inflaton in the first stage of the oscillating inflation for two limiting cases 
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2 is smaller than for 0g <C 2 because the potential is flat for 0g ^> 2 . As time 



Note that for 

passes and particles are produced, back reaction effects due to the increase of (5(j) 2 ) and (x 2 ) become relevant and the 
coherent oscillation of the tpo field is prevented by the change of the effective mass m^. 



In the frequency (3.5) of the Sifk field, since we consider the situation of < q < 1, the second term becomes 
negative in the region of 0q > 2 . Moreover, in the small vicinity of <pQ < 2 , this term rapidly grows and takes large 
positive value ~ M /4>c- These peculiar behaviors of the frequency lead to the enhancement of the 4> particle during 
the stage of the oscillating inflation. 

As for the x field, since the oscillation of the 4>o field is different from that of the ordinary inflaton potential, the 
increase of (x 2 ) occurs in a different way. In this case, although analytic investigations of the evolution of the x 
particle are rather difficult, we shall examine the parameter range of the critical value <p c and the coupling constant 
g where the x particle production takes place efficiently. 

With regard to the initial conditions of the fluctuation, we choose the states which correspond to the conformal 
vacuum as 



5<p k (0) : 

X k (0) 
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Sip k (0) = -iu Svk{0 - ) 6ip k (0), 
X k (0) = -icj Xh{o) X k (0), 



(3.9) 



(3.10) 



We investigate the evolution of the fluctuations of cj> and x fields with those initial conditions as the semiclassical 
problem. In order to clarify the situation we consider, we first investigate the case when the coupling g is zero (i.e. 
one-field case) and next investigate the case of g ^ (two-field case). 
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A. The case of g — 



In this subsection, we study the evolution of the inflaton quanta in the case of g = 0. This case was originally 
considered by Taruya [l9| including the metric perturbation. However, since his results are obtained neglecting the 
back reaction effect of created particles, we investigate how this effect would modify the growth of the fluctuation. 
The first thing we should notice is that the strength of the amplification of the fluctuation (5<p 2 ) strongly depends on 
the critical value <f> c . With the decrease of cf> Cl since the duration of the oscillating inflation becomes longer, parametric 
amplification of the <f> particle occurs more efficiently before the back reaction effect becomes significant. 

We investigate two concrete cases of <j> c = 10 _3 M p i and 4> c = 10~ 4 A/ p i with the value of q — 0.1 where the growth 
of the fluctuation can be expected. When q = 0.1, the initial value of <f>o is estimated as (j>o = 1.4 x 10 _2 M p i by 



Eq. ( [2.23D . 

Let us first consider the evolution of the 4>q field. In the case of <j> c = 10 _3 M p i, since the initial value of 4>o is larger 
than <j) c only by one order of magnitude, the number of e-foldings during the oscillating inflation is small [N = 0.93 



by Eq. ( 2. 21 )]. In Fig. 2, we depict the evolution of the 4>q field. The amplitude cf>o gradually decreases due to 
the expansion of the universe, and the oscillating inflation ends when 4>q is lowered to <p c . Numerical calculations 
indicate that this occurs at i = M 2 t/M p \ sa 0.37. After this, the <f>o field oscillates around the core region < (j) 2 , 
whose bare mass is given by = M 2 /cf> c . In the case of (f> c = 10 _4 M p i, since the oscillating inflation continues 
until <f>o w 10 M p i, we obtain the larger value of e-foldings N = 1.73 compared with the case of cj) c = 10 _3 M p i. 
Numerically, the oscillating inflation ceases at t « 0.45 in this case. 

With regard to the fluctuation of inflaton, the resonance band is very broad and the growth rate of the fluctuation 
becomes of order unity. We can expect the enormous amplification of the wave modes especially within the Hubble 
radius. We show in Fig. 3 the evolution of the 54>k field for two cases of k = k/(M 2 /M p \) = 0.1 and k = 100. We 
find that long-wave modes are not significantly enhanced compared with the modes inside the horizon scale. This 
result coincides with numerical calculations performed in Ref. JlSj ] which include the metric perturbation. As was 
presented in Ref. |19"| , the effect of the metric perturbation appears in the term 2k 2 (V/H) of the equation of the 
Mukhanov variable Q = 5<j) — lZ<fi/H, where 71 is the spatial curvature [See Eq. (8) in Ref. p9j]. However, since this 



term decreases faster than the term V^((/)o) in Eq. (2.9), we expect that adding this term to the equation of the 
fluctuation will not alter the evolution of the system. We have numerically checked that the growth of the fluctuation 
in Fig. 3 is almost the same as in the case where the 2k 2 (V/H) term is included. In the single-field case, there exists 
an exact solution for the Mukhanov variable in the long wavelength limit [See Eq. (9) in Ref. jL9|], and we can confirm 
that the amplitude of the fluctuation remains nearly constant. On the other hand, since the instability band is broad 
in the present model even for the large k, the momentum modes up to k 2 / a 2 < M 4 /(f> c 2 inside the Hubble horizon are 
effectively enhanced. In Fig. 3, we find that the fluctuation 5(f>k rapidly grows at the stage of the oscillating inflation 
for the case of k = 100. However, this increase is suppressed at i w 0.12, before the oscillating inflation terminates 
at i w 0.37. This is due to the fact that the back reaction effect of created particles becomes significant. We have 
numerically confirmed that the fluctuation continues to grow during the oscillating inflation if we neglect the back 
reaction effect as in Ref. |L9|. The increase of (5<fi 2 ) effectively changes both frequencies of ipo and Sipk fields as is 
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found by Eqs. (fTl) and (3.5), and the growth of the fluctuation finally stops when the (8(f> 2 ) terms become comparable 



to the preceding terms in Eqs. (3.4) and ( |3.5| ). It is expected that the back reaction effect becomes significant when 
(S(f> 2 ) increases of order <f> 2 . In the case of <f) c = 10 _3 M p i, numerical calculations show that the maximum value of 
the fluctuation is (8(f> 2 ) / « 10 _7 Mpj at i 0.12 (See Fig. 4). In this case, parametric amplification of the </> particle 
terminates when (6(j) 2 } increases up to (S(f> 2 ) / « 0.10;?. As is seen in Fig. 2, the coherent oscillation of the <fio field is 
hardly broken after t « 0.12 where the variance of the <j) particle reaches the maximum value. This means that the 
growth of (S(f> 2 ) affects the evolution of the 5<t>k field more strongly than that of the 0o field. 

With the decrease of </> c , since the oscillating inflation continues longer, the growth rate of the fluctuation in the first 
stage becomes larger. However, the final variance (5<fi 2 ) / is suppressed as 4> c becomes smaller. We depict in Fig. 5 the 
evolution of the fluctuation (S<fi 2 ) in the case of <p c = 10- 4 M pl . In this case, although the oscillating inflation continues 
until i w 0.4, the growth of the fluctuation stops much earlier: i w 0.05. Although the initial growth rate is larger 
than in the case of (j> c = 10 _3 M p i, the production of the </> particle itself terminates the growth of the fluctuation. 
This tendency is stronger with the decrease of </> c , and the maximum fluctuation becomes smaller. We have found the 
following relation for the typical value of <p c < 10 _3 M p i as 

(^)/«04 (3.11) 
Although the <f> particle production is possible at the initial stage of the oscillating inflation, the final variance is 



strongly suppressed with the decrease of </> c as Eq. (3.11). After the amplitude of the infiaton field drops under 



00 < 4>c, the system enters the ordinary reheating stage where the universe expands decelcratcdly. In this stage, the 
growth of the fluctuation (5(f> 2 ) can be no longer expected because the infiaton field behaves as the massive infiaton 
whose mass is = M 2 /cj) 2 . 

In the next subsection, we consider another field x coupled to infiaton and analyze how \ particles are produced 
during the oscillating inflation and subsequent oscillating phase by parametric resonance. 



B. The case of g 



Next, we consider the production of the \ particle coupled to the infiaton field. In the ordinary picture of the slow- 
roll inflation, the x particle production is inefficient during the inflationary phase. In the present model, however, 
since the infiaton field oscillates rapidly during the oscillating inflation, x particles can be produced in this stage by 
parametric resonance. 

The strength of resonance depends on the coupling g, the amplitude of c/)q (= <f>o), and the mass of the infiaton 
field (= m$). In the model of the scalar field x coupled to the massive infiaton <f>: V((f>, x) = m 2 (f> 2 /2 + g 2 cf> 2 x 2 /2, 
the equation of the x ne ld is reduced to the Mathieu equation [^6| at the linear stage. In this model, the resonance 
parameter q = g 2 <p 2 / (4m 2 ) is an important factor in determining whether the x particle production is efficient 
or not ||. When q is sufficiently large initially as qi 3> 1, parametric resonance turns on from broad resonance 
regimes. Although q decreases with the decrease of 4>o due to the cosmic expansion, the fluctuation (x 2 ) grows 
quasiexponentially until q drops under unity or the back reaction effect of created particles becomes significant. 
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In the oscillating inflation model, since the mass = + l^j continually changes during one oscillation 

of inflaton, the ipo held does not oscillate sinusoidally and the method based on the Mathieu equation is not valid 



except when 4>q drops under <p c . In the frequency (3.6) of the Xk field, since the last term is negligible when the x 
particle production occurs sufficiently, the equation of the \ field for the case of (S(f> 2 ) <C <Pq is approximately written 
by 



d " x 
dz^ Xk 



~ k ' 2 , ^2 , 9 2 <g 



rn 



i 2 x m 2 . 



X k « 0, (3.12) 



where z = m^t, k = fc/m^, and m x = m x /m ( p. When the initial value of g (po/ m 4, i s large, parametric resonance 



works out more efficiently. The initial value of </>o estimated by Eq. (2.23) for the typical value of q is smaller than 
the initial value 0o ~ 0.2-0.3M p i in the massive inflaton model. Moreover, as is found by Eq. (2.20), the amplitude of 
the 4>o field decays faster than in the model of the massive inflaton with tfio cx t . 



The effective mass of the 0o field is approximately written in two asymptotic regions as Eq. (3.8), and depends on 
the critical value <p c . The mass m$ ~ M 2 /4> 2 around <\>\ < cf> 2 plays the more important role than that of <\>\ > 4> 2 
for the development of the fluctuation, because \ particles are mainly produced in the vicinity of 4>q — where the 
0o field changes nonadiabatically. During one oscillation of the </>o field, gradually gets larger as (f>o approaches 
the minimum of the potential for the fixed value of <p c . Moreover, for the typical value of 4> c < 10 _3 M p i where the 
oscillating inflation occurs, = M 2 /(j> c is greater than the mass ~ 10 _6 M p i in the model of the massive inflaton. 
For example, in the case of (j) c — 10^ 3 M p i and q = 0.1, M 2 /4> c — 1.57 x 10~ 4 M p i; for <j) c = 10 _4 A/ p i and q = 0.1, 
M 2 /cj) c = 1.40 x 10~ 3 Af pl . 

In addition to the fact that cf>o decreases faster than in the model of the massive inflaton, the larger mass results 
in the restriction of the coupling constant g for an efficient particle production. Since the dependence of M 2 for 
the value 4> c is weak, the mass M 2 /<j) c monotonically increases with the decrease of <p c . This means that parametric 
resonance does not take place unless we choose a large coupling constant g as <j> c decreases. When <f) c is fairly large as 
(j> c > 0.1M p i, M 2 /<fi c becomes comparable to the mass of the massive inflaton model, but in this case the creation of 
4> particles can not be expected because of the absence of the oscillating inflationary phase. What we are interested 
in is the case of 4> c < 10 _3 M p i where the oscillating inflation occurs and both of <j> and \ particles arc generated. In 
what follows, we examine the x particle production in two cases of <fi c — 10~ 3 AT p i and <f> c = 10~ 4 Af p i for the massless 
X particle, and add the discussion of the case where the mass of the x particle is included at the final of this section. 

When <fi c — 10 Afpi, numerical calculations indicate that the increase of (x 2 ) can take place for g > 5 x 10~ 3 . 
However, parametric resonance is weak for the case of g < 0.03. As compared with the model of V(<p, x) = rn 2 (j) 2 /2 + 
g 2 (fi 2 X 2 /2 in which the x particle production takes place for g > 10~ 4 ||, larger values of g are required even in the 
case of 4> c = 10 _3 M p i where the oscillating inflation marginally occurs. Since the fluctuation (x 2 ) does not grow well 
for g < 0.03, the evolution of the </>o field is hardly affected by the x particle production. For the case of g > 0.05, 
the back reaction effect of the created x particle begins to be significant. We depict in Fig. 6 the evolution of (x 2 ) 
for two cases of g = 0.03 and g = 0.07. Although parametric resonance evidently occurs for g — 0.03, this is rather 
inefficient and (x 2 ) takes the maximum value (x 2 ) f = 2.3 x 10 _11 A/ 2 j at t = 0.177. In this case, the x field deviates 
from the resonance bands with the decrease of 0o due to the cosmic expansion before the oscillating inflation ceases 
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at i w 0.37. Since the maximum fluctuation of (54> 2 ) is almost the same as the g = case, the back reaction effect of 
4> particles onto the 4>o field can be marginally negligible. Hence in the case of g < 0.03, the coherent oscillation of the 
4>o field is hardly prevented by the back reaction effects of both <\> and x particles. When g — 0.07, the evolution of 
the system shows different characteristics compared with the case of g < 0.03. In Fig. 7, the evolution of the (j>o field 
for the g = 0.07 case is depicted. We find that the coherent oscillation of the <fio field is broken for i > 0.12, which 
is different from the 3 = case in Fig. 2. As for the <j> particle, the maximum fluctuation is (S(j) 2 ) / = 2.7 x 10 _6 M 2 [ 
at t — 0.12. This maximum variance is greater than in the case of g = by one order of magnitude (See Fig. 8). 
As was noticed in the previous subsection, the growth of (6(f) 2 ) stops when the contribution of the third term in 
the frequency (3.5) becomes comparable to the second term. Since the third term mostly takes negative values, the 
increase of the g 2 (x 2 ) term due to the \ particle production generally assists the <ft particle production to continue 
longer. This results in the larger value of the final fluctuation (S<fi 2 ), and the back reaction effect due to the <j> particle 
production becomes relevant at i m 0.12. This behavior is found in Fig. 7. However, since the 4>q field still oscillates 
after the back reaction effect of the <j) particle becomes important, the x particle continues to be enhanced after 
(« 0.12. The fluctuation of (x 2 ) does not increase significantly in the initial stage of the oscillating inflation, because 
the expansion rate of the universe is large, (x 2 ) begins to increase rapidly after t ~ 0.15, and reaches the maximum 
value (x 2 ) f — 2.4 x 10~ 7 M 2 j at i = 0.28. The back reaction effect of the x particles as well as the decrease of <fro 
terminates the parametric amplification of x particles. When g = 0.1, the increase of (x 2 ) is more significant and 
the final variance is (x 2 ) f = 3.3 x 10~ 6 Afp[ . In this case, since the maximum variance is larger than in the case of 
g = 0.07 by one order of magnitude, the coherent oscillation of the </>o field is more strongly prevented by the growth 
of (x 2 )- When g > 0.3, the final variance is suppressed because the back reaction effect becomes quite significant. We 
find that (x 2 )/ takes the maximum value (x 2 )max ~ 10~ 5 M 2 ! for g w 0.3 (See Fig. 9). 

As for the momentum modes of the produced x particle, we show the evolution of the Xk field in two cases of k = 
and k — 100 for g — 0.1 in Fig. 10. Since the larger k makes the Xk field deviate from the resonance band in Eq. ( p,12j ), 
higher momentum modes of the x particle are not sufficiently produced. The low momentum modes mainly contribute 
to the growth of the fluctuation (x 2 )- This property is the same as in the model of V(<j), x) = m 2 (j) 2 /2 + g 2 <p 2 x 2 /2 m- 
As was found in Ref. |3(if| , there is a possibility that long wavelength modes of the metric perturbation are amplified 
considering the perturbed metric. This issue will be discussed elsewhere j27j. 

In the case of <fi c — 10 M p i, we need further large values of g to yield an efficient resonance. Numerically, we 
find that the coupling is required g > 0.01 for the development of the x fluctuation. Even when g = 0.1, parametric 
resonance is not so efficient compared with the case of <j> c — 10 _3 M p i. We depict in Fig. 11 the evolution of (x 2 ) for 
g = 0.07, 0.1, 0.5 cases respectively. When g — 0.07, the maximum value of the fluctuation is (x 2 ) / = 1.2 x 10 _10 Af 2 



which is much smaller than in the case of <fi c = 10 3 M p \ with the same value of g. As the coupling increases further, 



V' 

V " ' ' ' '' " " " 

we obtain larger values of the maximum fluctuation (x 2 )/- However, for g > 0.5, the back reaction effect of created 

X particles becomes significant and the final variance is suppressed. In the case of <f> c — 10 _4 M p i and g — 0, the 

fluctuation of the field rapidly grows in the initial stage and reaches the maximum value (5<fi 2 ) f w lO" 9 ^^ at 

i = 0.05. Even if we take into account the interaction with the x field, this maximum variance is hardly altered 

because the increase of (x 2 ) is weak for t < 0.05 as is found in Fig. 11. The back reaction effect of <f> particles onto 
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the 4>q field is marginally negligible as the case of g — 0. When g = 0.1, the final variance is (x 2 ) / = 2.5 x 10~ 9 M 2 l5 
and the growth of the fluctuation hardly affects the evolution of the </>o field. However, in the case of g = 0.5, (x 2 )/ 
increases up to (x 2 ) / = 1-3 x 10~ 6 M 2 j at i w 0.17. In Fig. 12, we find that the increase of x particles prevents the 
coherent oscillation of the 4>q field for i > 0.17. When g > 0.5, since the \ particle production stops by the back 
reaction effect, the final variance gradually decreases with the increase of g. As a result, the maximum variance for 
the case of (j) c = 10~ 4 M pl is (x 2 ) max ~ IQ'^M^ for g « 0.5. 

With the decrease of <f> c ( < 10 _5 M p i), the development of the \ fluctuation is not expected unless g is unnaturally 
large. Moreover, in this case, although the 4> particle production is possible in the initial stage of the oscillating 
inflation, the final fluctuation (8(f) 2 ) / is strongly suppressed. This means that it is difficult to obtain the sufficient 
amount of <fi and x particles with the natural coupling g < 1 in the case of 4> c < 10 _5 M p i. For large critical values 
4> c > 10~ 3 A/ p i, the x particle production occurs effectively for the natural coupling. With the increase of </> c , the 
structure of resonance for the x field approaches that of the model V(4>, x) = m 2 (f) 2 /2+g 2 (f> 2 x 2 and the investigation 
based on the Mathieu equation performed in Ref. j^] can be applied. Since the mass around the minimum of the 
inflaton potential becomes smaller, the x particle production occurs efficiently. In this case, however, since the duration 
of the oscillating is too short, the <j) particle production is hardly expected. In the case of the massless x particle, we 
conclude that both </> and x particles are most efficiently created in the parameter range <fi c = 10 _3 -10~ 4 Mpi for the 
natural coupling g. 

Finally, we discuss the case where the mass of the x particle is included. We expect that the large mass makes the 
X field deviate from the resonance band by the relation (3.12). Let us consider the case of <fr c = 10 _3 M p i and q = 0.1. 
In the massless case, the x particle production occurs relevantly for g > 0.05. If the mass of the x particle is taken 
into account, numerical calculations show that heavy particles up to the order of m x = 10 14 GeV can be produced 
enough for g > 0.05. In Fig. 13, the evolution of the fluctuation (x 2 ) is depicted in two cases of m x — 3 x 10 14 GeV 
and m x — 1 x 10 15 GeV for g — 0.05. The final variance (x 2 ) / ~ lO -8 ^/^ with mass m x — 3 x 10 14 GeV is almost the 
same as in the case of the massless x particle with the same coupling g. This means that the massive x particle whose 
mass is m x < 3 x 10 14 GeV is produced with almost the same amount as the massless x particle for g = 0.05. However, 
the achieved amount of the produced x particle decreases with the increase of m x for m x >3x 10 14 GeV. In the case 
of g = 0.05, although the x particle with mass m x < 7 x 10 14 GeV can be created, it is difficult to generate the massive 
particle whose mass is m x > 1 x 10 15 GeV (See Fig. 13). In the model of V(<p, x) = m 2 2 /2 + g 2 2 x 2 /2, analytic 
investigations show that the x particle which satisfies the condition m x < ^fg 10 15 GeV can be generated p,0. When 
m x — 10 14 GeV, numerical calculations performed in Ref. || reveal that the x particle production occurs for g > 0.06 
p8| . This condition is almost the same as in the oscillating inflation model in spite of the difficulty in producing 
the massless x particle without the larger g compared with the V((f>, x) = m 2 </> 2 /2 + <? 2 </> 2 x 2 /2 model. The reason is 
simple. In Eq. fl3.12| ), the massive x particle production occurs when the initial value of the resonance term g 2 <^lmh 
is much greater than unity [(g 2 4>o/'n r i'i)i ^ 1] and m 2 is much smaller than this term [to 2 <C {g 2 <t>a/'rn 2 c j > )i\- In the case 
of <f> c = 10" 3 M p i, q = 0.1, and g = 0.05, since M = 3.96 x 10~ 4 M pl by Eq. ( |2.2fl ), we obtain {g 2 <Pl/m 2 > ) i w 4 x 10 3 . 
Although 4>q decreases by the cosmic expansion and m$ becomes large in the vicinity of 4>o = 0, the massive x particle 
production is possible as long as ff 2 ^o/ TO l ^ 1 an d ^ ff 2 ^o/ m I- Since the mass in this case is estimated as 
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TTigp w 10 _5 M pl for 4>l > 4> 2 C and w 10~ 4 M pl for 4>o < <^c by Eq. this is heavier than the mass m w 10~ 6 M p i 

in the U(</>, x) = m 2 </> 2 /2 + g 2 <fr 2 x 2 /2 model. Even in the case where the normalized mass is m x = 1, this corresponds 
to the rather heavy particle whose mass is m x > 10 14 GeV. For the typical value <j> c < 10 _3 M p i in the oscillating 
inflation model, since is at least by one order larger than m = 10 Af p i, the massive x particle whose mass is of 
order m x = 10 GeV can be generated for the coupling g where the massless x production occurs. As the coupling 
g increases, more massive particles are generated. In the case of 4> c = 10 _3 M p i, we find that the production of the 
massive particle whose mass is of order m x — 10 15 GeV relevantly occurs for the coupling g > 0.2. Even the GUT 
scale bosons m x = 10 16 GeV are produced in the initial stage of the oscillating inflation for the case of g = 1, although 
the achieved amount is small (See Fig. 14). 

With the decrease of <j> c , the larger values of g are required in order to produce the massive x particle as is 
the same with the massless x particle. However, in the parameter range of g where the massless x particles are 
sufficiently produced, we can also expect the generation of the superheavy particles needed for the success of the GUT 
baryogenesis. For example, in the case of <p c = 10~ M p \, numerical calculations indicate that x particles up to the 
mass m x < 10 14 GeV and m x < 10 15 GeV are created for the coupling g > 0.07 and g > 0.3 respectively. It will be 
interesting to investigate how these superheavy bosons produced in the oscillating inflation model would affect the 
scenario of the GUT baryogenesis. 



IV. CONCLUDING REMARKS AND DISCUSSIONS 



In this paper we have considered parametric amplification of a scalar field x coupled to an inflaton field <j) as well 
as the enhancement of the fluctuation of the <p field in the oscillating inflation model. Since the inflaton potential 
V(4>) is the non-convex type where d 2 V/d(f> 2 takes negative values in the regions not too far from the minimum of 
the potential, inflation is realized even during the oscillation of the inflaton field. The oscillating inflation model is 
characterized by the critical value of <j> c which indicates the scale of the core part of the potential. Inflation takes 
place while the <\> field moves in the regions of \<$>\ > 4> c . With the decrease of <j) c , the total amount of inflation 
becomes larger, because the duration during which the oscillating inflation takes place becomes longer. However, it 
is rather difficult to obtain the sufficient number of e-foldings such as TV > 60 only by the stage of the oscillating 
inflation, and the main contribution to the number of e-foldings is achieved by the ordinary slow-roll inflation which 
precedes the oscillating inflation. We have estimated the energy scale of the oscillating inflation M 4 by fitting the 
primordial density perturbation observed by COBE. The oscillating inflation terminates when the amplitude of the 
4>q field decreases to of the order <f> c by the expansion of the universe. 

In this model, the fluctuation of the (f> field grows during the oscillating inflation due to the existence of the non- 
convex part of the potential. We examined the growth of the fluctuation for two typical cases (j> c — 10 -3 M p i and 
4> c = 10~ 4 M p i including the back reaction effect of created particles based on the Hartree approximation. We have 
found that the wave modes inside the Hubble horizon are effectively enhanced, while the growth of long-wave modes is 
small as was shown in Ref. [pi. Although the growth rate of the fluctuation (84> 2 ) in the initial stage of the oscillating 
inflation becomes larger with the decrease of <j> c , the final fluctuation is significantly suppressed. This indicates that 



13 



the generation of the <\> particle itself plays a crucial role for the termination of parametric resonance. As for the final 
variance of the <j> particle, we find the relation (6<p 2 ) / ss O.ltfi 2 for the case of <fi c < 10 _3 M p i. 

With regard to the x held coupled to inflaton, the enhancement of the fluctuation mainly occurs when the <f>Q field 
moves nonadiabatically around the minimum of the potential. Since the mass for |0q| < 4> c is approximately 
written as w M 2 /0 C , it is larger than that of the massive inflaton of the chaotic inflation model for the value of 
(j) c < 10 _3 M p i. As m<t, becomes larger, this restricts the effective x particle production. Moreover, the amplitude of the 
homogeneous inflaton field decreases faster than in the model of the massive inflaton. This means that the x particle 
production is not possible unless we select rather large values of the coupling constant g. This tendency becomes 
stronger with the decrease of <j> c . For example, in the case of 4> c — 10 _3 Af p i and q = 0.1, the coupling g > 0.05 is 
required to yield the sufficient x particle production, the value of which is two orders of magnitude larger than in the 
model of V(<j>, x) = m 2 <p 2 /2 + g 2 <p 2 x 2 /2. For g > 0.05, the back reaction effect onto the </>o field due to the x particle 
production is significant, and finally terminates parametric resonance. The final variance is suppressed by the back 
reaction effect with the increase of g, and the maximal variance for the <fi c = 10~ 3 M p i case is (x 2 ) / ~ 1Q~ 5 M 2 ^ when 
g « 0.3. In the case of 4> c = 10~ 4 M p i, we need g > 0.1 for an effective resonance, although the x particle production 
occurs for g > 0.01. The maximal variance for the 4> c = 10~ 4 M p i case is (x 2 ) f ~ 10~ 6 A/pj when g w 0.5. As </> c 
decreases further, we no longer expect the growth of the x fluctuation with the natural coupling g < 1. Moreover, 
since the final fluctuation of the 4> particle is also suppressed, the small 4> c such as cf> c < 10~ 5 M p i is not favorable for 
the production of a sufficient amount of <p and x particles. 

In the oscillating inflation model, the superheavy x particle with mass greater than m x = 10 14 GeV can be copiously 
produced. For example, in the case of <\> c = 10~ 3 M p i, we find that the x particle whose mass is of order 10 14 GeV 
and 10 15 GeV is effectively enhanced for g > 0.05 and g > 0.2 respectively. The GUT scale boson m x ~ 10 16 GeV 
can be also generated in the initial stage if the coupling g is of order unity, although the final amount is small. In 
Ref. |22j , the authors considered the GUT scale baryogenesis with the massive inflaton assuming that some of the 
initial inflaton energy is efficiently transferred to the bosons with the mass 10 14 GeV. These massive bosons decay 
into lighter particles, and produce the net baryon number. Although we do not consider such a decaying process in 
this paper, it is of interest how the production of the superheavy particle would affect the baryon asymmetry in the 
universe. 

In this paper we make use of the Hartree approximation, which is essentially the mean field approximation. This does 
not include the rescattering effect which becomes important as x particles are sufficiently produced. The scattering 
between <fi and x particles may reduce the final amount of fluctuations. For a complete study including the nonlinear 
effect of the particle production, we should perform the lattice simulations and compare them with the mean field 
approximation performed in this paper. 

Finally, we comment on the case where the metric perturbation is included. In the single field case, even if we 



consider the effect of the metric perturbation, the V^(<f>o) term in Eq. (2.9) is much more important than the 
2k 2 (V/H) term which is the gravitational origin. Hence the situation is almost the same as the case when the metric 
perturbation is neglected. As we have showed, the long-wave modes are not significantly enhanced compared with 
the modes inside the Hubble Horizon in the single-field case. This result is consistent with other inflation models 
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in the single-field case p9fl . In the two-field case, it was suggested that super-Hubble metric perturbations can be 
amplified in broad classes of models p0| . In the present model, since the \ fluctuation of super-Hubble modes will 
be exponentially suppressed for large values of g > 1 |sT| , we may not expect the enhancement of super-Hubble 
metric perturbations significantly. However, it is worth investigating to investigate these issues including mode-mode 
coupling between field and metric perturbations, because this may lead to some imprints on the spectrum of density 
perturbations. These issues are under consideration. 

ACKOWLEDGEMENTS 

The author would like to thank Bruce A. Bassett, Kei-ichi Maeda, Atsushi Taruya, Takashi Torii, and Hiroki 
Yajima for useful discussions. This work was supported partially by a Grant-in- Aid for Scientific Research Fund of 
the Ministry of Education, Science and Culture (No. 09410217 and Specially Promoted Research No. 08102010), and 
by the Waseda University Grant for Special Research Projects. 



[1] A. H. Guth, Phys. Rev. D 23, 347 (1981); K. Sato, Mon. Not. R. Astron. Soc. 195, 467 (1981); Phys. Lett. 99B, 66 (1981). 
[2] E. W. Kolb and M. S. Turner, The Early Universe (Addison- Wesley, Redwood City, California, 1990); A. D. Linde, Particle 

Physics and Inflationary Cosmology (Harwood, Chur, Switzerland, 1990). 
[3] A. D. Dolgov and A. D. Linde, Phys. Lett. 116B, 329 (1982); L. F. Abbott, E.Fahri, and M. Wise, Phys. Lett. 117B, 29 

(1982). 

[4] J. Traschen and R. H. Brandenberger, Phys. Rev. D 42, 2491 (1990); Y. Shatanov, J. Trashen, and R. H. Brandenberger, 

Phys. Rev. D 51, 5438 (1995). 
[5] L. Kofman, A. Linde, and A. A. Starobinsky, Phys. Rev. Lett. 73, 3195 (1994). 

[6] S. Khlebnikov and I. I. Tkachev, Phys. Lett. B390, 80 (1997); S. Khlebnikov and I. I. Tkachev, Phys. Rev. Lett. 79, 1607 
(1997). 

[7] T. Prokopec and T. G. Roos, Phys. Rev. D 55, 3768 (1997); B. R. Greene, T. Prokopec, and T. G. Roos, Phys. Rev. D 
58, 6484 (1997). 

[8] L. Kofman, A. Linde, and A. A. Starobinsky, Phys. Rev. D 73, 3258 (1997). 

[9] D. Boyanovsky, H. J. de Vega, R. Holman, D. S. Lee, and A. Singh, Phys. Rev. D 51, 4419 (1995); D. Boyanovsky, M. D' 
Attanasio, H. J. de Vega, R. Holman, D. S. Lee, and A. Singh, Phys. Rev. D 52, 6805 (1995); D. Boyanovsky, H. J. de 
Vega, R. Holman, D. S. Lee, and A. Singh, J. F. J. Salgado, Phys. Rev. D 54, 7570 (1996); D. Boyanovsky, D. Cormier, H. 
J. de Vega, R. Holman, A. Singh, and M. Srednicki, Phys. Rev. D 56, 1939 (1997); D. Boyanovsky, R. Holman, S. Prem 
Kumar, Phys. Rev. D 56, 1958 (1997). 

[10] D. I. Kaiser, Phys. Rev. D 53, 1776 (1995); ibid D 56, 706 (1997). 

[11] S. Khlebnikov and I. I. Tkachev, Phys. Rev. Lett. 77, 219 (1996). 

[12] D. T. Son, Phys. Rev. D 54, 3745 (1996). 

[13] J. Baacke, K. H eitmann, C. Patzold, Phys. Rev. D 55, 2320 (1997); ibid D 56, 6556 (1997); J. Baacke and C. Patzold, 



hep-ph/9906417. 



[14] P. B. Greene, L. Kofman, A. Linde, and A. A. Starobinsky, Phys. Rev. D 56, 6175 (1997). 
[15] J. Garcia-Bellido and A. Linde, Phys. Rev. D 57, 6075 (1998). 

[16] S. Tsujikawa, K. Maeda, and T. Torii, Phys. Rev. D 60, 123505 (1999); see also B. A. Bassett and S. Liberati, Phys. Rev. 
D 58, 02130 2 (1998); S. Tsuj ikawa, K. Maeda, and T. Torii, Phys. Rev. D 60, 063515 (1999); S. Tsujikawa, K. Maeda, 



and T. Torii, 



hep-ph/9910214| , to appear in Phys. Rev. D. 



[17] T. Damour and V. F. Mukhanov, Phys. Rev. Lett. 80, 3440 (1998). 
[18] A. R. Liddle and A. Mazumdar, Phys. Rev. D 58, 083508 (1998). 
[19] A. Taruya, Phys. Rev. D 59, 103505 (1999). 

[20] L. Kofman, A. Linde, and A. A. Starobinsky, Phys. Rev. Lett. 76, 1011 (1996); I. I. Tkachev, Phys. Lett. B376, 35 (1996); 
S. Kasuya and M. Kawasaki, Phys. Rev. D 56, 7597 (1997). 



15 



[21] S. Khlebnikov, L. Kofman, A. Linde, and I. I. Tkachev, Phys. Rev. Lett. 81, 2012 (1998); I. I. Tkachev, S. Khlebnikov, L. 

Kofman, A. Linde, Phys. Lett. B440, 262 (1998); S. Kasuya and M. Kawasaki, Phys. Rev. D 58, 083516 (1998). 
[22] E. W. Kolb, A. D. Linde, and A. Riotto, Phys. Rev. Lett. 77, 3716 (1996); E. W. Kolb, A. Riotto, and I. I. Tkachev, Phys. 

Lett. B423, 348 (1998). 

[23] In this paper we do not consider the nonminimal coupling of the or x field with the scalar curvature R. However, there 

is a possibility that the nonminimal coupling would lead to the gravitational particle production as in Ref. Jlq ]. 
[24] D. Cormier and R. Holman, Phys. Rev. D 60, 041301 (1999). 

[25] D. H. Lyth, Phys. Rev. D 31, 1792 (1985); A. R. Liddle and D. H. Lyth, Phys. Rep. 231, 1 (1993). 
[26] N. W. Mac Lachlan, Theory and Applications of Mathieu Functions (Dover, New York, 1961). 
[27] S. Tsujikawa, in preparation. 

[28] This criterion was also obtained by analytic estimates performed in Ref. j^]. 

[29] H. Kodama and T. Hamazaki, Prog. Theor. Phys. 96, 949 (1996); Y. Nambu and A. Taruya, Prog. Theor. Phys. 97, 83 
(1997); F. Finelli and R. Brandenberger, Phys. Rev. Lett. 82 1 362 (1999); M. Parry and R. Easther, Phys. Rev. D 59 
061301 (1999); S. Tsujikawa and B. A. Bassett, tiep-ph/0003068 . 

[30] B. A. Bassett, D. I. Kaiser, and R. Maartens, Phys. Lett. B455, 84 (1999); B A. Bassett, F. T amburini, D. I. Kaiser, and 



R. Maartens, Nucl. Phys. B561, 188 (1999); B. A. Bassett and F. Viniegra, |hep-ph/9909353 



B. A. Bassett, C. Gordon, 



R. Maartens, and D. I. Kaiser, Phys. Rev. D. 061302(R) (2000); F. Finelli and R. Brandenberger, ^iep-ph/0003172 
[31] P. Ivanov, Phys, Rev. D 61, 023505 (2000); K. Jedamzik and G. Sigl, Phys, Rev. D 61, 023519 (2000). 



1G 



Figure Captions 



FIG. 1: The potential for q = 0.1 in the oscillating inflation 

model, which is composed of the non-convex region 
\(f)\ > 4> c and the core region \<f>\ < <fi c . 

FIG. 2: The evolution of the <j>o field in the case of <p c — 

10- 3 M p i, q = 0.1, and g = 0. The stage of the 
oscillating inflation starts from ^ w 1.4 x 10~ 2 M p i, 
and ends at i w 0.37 when the amplitude of the </>o 
field becomes of order <p c . 

FIG. 3: The evolution of the real part of the fluctuation 5(j)k 

for two momentum modes of k = 0.1 and k = 100 
in the case of 4> c = 10 _3 M p i, q = 0.1, and g = 0. 
For the long wave mode k — 0.1, the growth of the 
fluctuation is weak, but for the mode k = 100, the 
fluctuation grows rapidly until the back reaction ef- 
fect becomes significant at t w 0.12. 

FIG. 4: The evolution of (S(f> 2 ) in the case of <\> c = 10~ 3 M p i, 

q = 0.1, and g = 0. The fluctuation reaches the 
maximum value (S(f> 2 ) f w lO^M^ at t w 0.12. 

FIG. 5: The evolution of (5(f) 2 ) in the case of 4> c = 10~ 4 M pl , 

q = 0.1, and g = 0. Although the initial growth 
rate is larger than in the case of <j> c = 10~ 3 M p i, 
the fluctuation soon reaches the maximum value 
(5<j) 2 ) f w 8.0 x 10- 10 M^, which is by two orders 
smaller than in the case of <j> c = 10 _3 M p i. 

FIG. 6: The evolution of the fluctuation (x 2 ) in the case of 

(j) c = 10- 3 M pl and q = 0.1 for g = 0.03 (bottom) and 
g = 0.07 (top). For g = 0.03, the maximum value 
of the fluctuation is (x 2 ) / = 2.3 x 10" n M 2 j; and for 
g = 0.07, (x 2 )/ - 2.4 x 10"^. 

FIG. 7: The evolution of the </>o field in the case of <j> c = 

10 _3 M p i, q = 0.1, and g = 0.07. The coherent oscil- 
lation is a bit broken at i w 0.12 due to the increase 
of (S(f) 2 }. With the increase of (x 2 ), this also changes 
the frequency of the </>o field at i w 0.28. 

FIG. 8: The evolution of (6(f> 2 ) in the case of (f> c = 10" 3 M pl , 

q = 0.1, and g = 0.07. With the growth of (x 2 ), 
this assists the <f> particle production, and the final 
variance (5<j) 2 ) / = 2.7 x 10~ 6 M 2 [ becomes larger than 
in the case of <p c = 10 _3 M p i and g — by one order 
of magnitude. 

FIG. 9: The final variance (x 2 )/ as the function of g in the 

case of <p c = 10~ 3 M p i and q = 0.1. We find that 
(x 2 )/ takes the maximum value (x 2 )max ~ 10 _5 M 2 1 
when g w 0.3. 
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FIG. 10: The evolution of the real part of the fluctuation Xk 
for two momentum modes of k — and k — 100 in 
the case of <j) c =_10 _3 M p i, q = 0.1, and g = 0.1. The 
long wave mode k = is strongly enhanced compared 
with the wave mode k = 100. 

FIG. 11: The evolution of the fluctuation (x 2 ) in the case of 
<j> c = 10" 4 Mpi, q = 0.1, and for g = 0.07 (bottom), 
g = 0.1 (middle), and g — 0.5 (top) respectively. 
The final variances are (x 2 ) / = 1.2 x 10- 10 M 2 1; 2.5 x 
lO^M^, and 1.3 x lO^M^ respectively. 

FIG. 12: The evolution of the (j>o field in the case of <\> c = 
10" 4 Af p i, q = 0.1, and g = 0.5. While the back 
reaction effect of <p particles can be negligible, the 
increase of (x 2 ) strongly affects the evolution of the 
00 held from t w 0.17. 

FIG. 13: The evolution of (x 2 ) for <j) c = 10~ 3 Mpi, q = 0.1, and 
g = 0.05 in two cases of m x = 3x 10 14 GeV (top) and 
m x = 1 x 10 15 GcV (bottom). The x particle whose 
mass is of order m x — 10 14 GeV is generated, but the 
production of the massive particle with m x > 10 15 
GeV is hardly expected for g = 0.05. 

FIG. 14: The evolution of (x 2 ) for C = lO^Afpi and q = 0.1, 
in two cases of g = 0.5, m x = 1 x 10 15 GeV (top); 
and g = 1, m x = 1 x 10 16 GeV (bottom). In the 
case of g = 0.5, the \ particle whose mass is of order 
m x — 10 15 GeV is effectively enhanced. In the case 
of g — 1, the GUT scale boson m x — 10 16 GeV is a 
little produced in the initial stage of the oscillating 
inflation. 
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